Configurations and Observables in an Ising Model with Heat Flow 
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We study a two dimensional Ising model between thermostats at different temperatures. By ap- 
plying the recently introduced KQ dynamics, we show that the system reaches a steady state with 
coexisting phases transversal to the heat flow. The relevance of such complex states on thermo- 
dynamic or geometrical observables is investigated. In particular, we study energy, magnetization 
and metric properties of interfaces and clusters which, in principle, are sensitive to local features of 
configurations. With respect to equilibrium states, the presence of the heat flow amplifies the fluctu- 
ations of both thermodynamic and geometrical observables in a domain around the critical energy. 
The dependence of this phenomenon on various parameters (size, thermal gradient, interaction) is 
discussed also with reference to other possible diffusive models. 



PACS numbers: 05.60.-k Transport Processes - 05.50.+q Lattice Theory and statistics 
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I. INTRODUCTION 

The study of systems undergoing heat flows is a classi- 
cal topic in non equilibrium statistical mechanics. Several 
important results have been obtained, especially for one 
dimensional models with continuous symmetries, such as 
chains of anharmonic oscillators (see e.g. [l[ for a review). 
On the contrary, there are very few results for discrete 
models in two dimensions. A ferromagnetic rectangular 
Ising lattice with a "cylindrical" geometry, i.e. opposite 
borders at temperatures T\ and T 2 in one direction, and 
periodic conditions in the other one, has been introduced 
in by Harris and Grant, and in [H by Saito, Take- 
sue and Miyashita (see also [J] for recent developments 
on related matter). However, severe restrictions on the 
admitted temperature intervals were present in both pa- 
pers, due to intrinsic limitations of the microcanonical 
dynamics used there (Creutz or Q2R rules). 

Such restrictions have been removed in Q by intro- 
ducing a peculiar new dynamics, briefly denoted as "KQ 
dynamics" , combining the advantages of the Q2R and 
Kadanoff-Swift rules. In this way, due to an effective 
ergodicity in the whole range of temperature intervals 
(Ti,T 2 ), steady states take place for all imposed tem- 
peratures. In particular, for T\ < T c < T 2 (where T c 
denotes the equilibrium critical temperature), different 
phases steadily coexist: a magnetized phase near the cold 
border at T\, a paramagnetic phase near the opposite hot 
border at T 2 , and an intermediate phase around the re- 
gion at energy density E c , the mean energy correspond- 
ing, at equilibrium, to T c . Moreover, the transport prop- 
erties of the system are well described by introducing an 
energy dependent diffusivity. This occurs in a smooth 
way, possibly except around E c , where the specific heat 
diverges and the diffusivity vanishes in the thermody- 
namic limit. 

Transport apart, an open problem - and our main item 
indeed - is the physical relevance of such steady states, 



characterized by many coexisting phases, as they are dis- 
tinct from homogeneous equilibrium states. More pre- 
cisely, for local physical observables, we ask if a portion 
of the cylinder has recognizable and peculiar properties 
when a heat flow passes through it. In particular, we 
shall concentrate on sections perpendicular to the flow 
(columns, vertical bands). Two kinds of physical observ- 
ables will be considered: thermodynamic quantities, such 
as energy density and magnetization, and geometrical- 
dynamical observables, for which the role of the configu- 
rations driven by the dynamics is predominant. The lat- 
ter observables are based on the metric properties of the 
configurations, which may involve very different items: 
the integral of pointwise differences (i.e. the well known 
Hamming distance) , which in some cases assumes an "en- 
ergetic" meaning, or the measure of differences in clus- 
ter distributions (Rohlin distance), an information-based 
metrics requiring the formalism of partition spaces. 

The main point is the existence of an energy band AE, 
starting just below E c , where the observable fluctuations 
are remarkably wider for a system undergoing a heat flow 
with respect to thermalized or close systems. The same 
happens to the distances between configurations. All this 
may be read as evidence of a larger variability of the 
system when it is far from equilibrium. These features 
strongly depend on the size L, and they disappear in 
the thermodynamic limit L — > oo. More precisely, as 
expected, they vanish as soon as the energy gradient be- 
tween neighbouring columns becomes infinitesimal and 
local equilibrium is reached. However, since real systems 
are characterized by finite gradients and finite sizes, such 
large fluctuations could be relevant in the study of meso- 
scopic systems with stationary flows. 

We recall that there are examples of exotic dynamics 
where the local equilibrium is not reached even for in- 
finitesimal gradients [||. Remarkably, also in such cases 
fluctuations are larger in the presence of heat flow. 

A number of questions arise. For instance, how much 
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do these features depend on the chosen dynamics? And 
which is the role of the specific spin interaction? As for 
the former question, the robustness of our results has 
been tested by many checks, improving, in addition, the 
reliability of the results described in [5[ . The latter ques- 
tion is evidently crucial for the possible physical relevance 
of the results. Now, for a purely diffusive process, e.g. a 
Random Walk (RW), analogous experiments clearly in- 
dicate the absence of the described phenomena, showing 
the essential role played by the interaction. However, a 
deeper insight on the nature of admissible interactions 
would require a more sophisticated analysis, not devel- 
oped here. The same holds for the role of other possible 
relevant parameters, such as the topology of the underly- 
ing structure or the presence of noise in the interactions. 

The paper is organized as follows: in Section [IT] the 
model is introduced, with KQ dynamics (jll Ap . and with 
definitions and notations for the quantities involved in ex- 
periments plB llrTC]) ; in Section Hill we review the main 
results obtained from numerical experiments. Problems 
recalled above (relevance of KQ dynamics on the results, 
etc.) are discussed in Section HV1 with further comments 
and perspectives on future work. Finally, in the Ap- 
pendix, we summarize the essential information on the 
formalism necessary to define the Rohlin distance in par- 
titions spaces. 



II. MODEL, DYNAMICS, NOTATIONS 

A. The Cylindrical Ising Model 

The cylindrical Ising model considered in [3[ and Q is 
a Lx x Ly rectangular lattice, with periodic conditions in 
the Y direction and open boundaries in the X direction. 
We assume Lx = Ly = L. The spin variable <j XtV may 
be 1 or —1, and adjacent opposite spins give an energy 
unit to the system. Thus, by denoting {x, y) the nearest 
neighbours of (x,y), the normalized total energy E to t is: 

P 1 \ " \ " 1 <J x,yO'{x,y) , . 

A * ot " Itfl^l^ 2 ' [ ] 

x -v (x,y) 

The lattice is naturally sliced into "columns" with a 
circular symmetry. The first and last columns, i.e. the 
left and right borders, interact with two thermostats, 
simulated by two sets of supplementary columns evolving 
with the usual equilibrium Metropolis algorithm (see [5[ 
for details). The Boltzmann's constant K is assumed to 
be 1. 

Internal sites must evolve preserving the energy, and 
the microcanonical rule used throughout the paper is the 
KQ dynamics introduced in 0] , for the reasons discussed 
there. In order to define such a dynamics, we must pre- 
viously recall the Q2R and Kadanoff-Swift (KS) moves: 
Q2R move: in every chosen site the spin is forced to 
flip whenever energy is preserved, i.e. when half spins 
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FIG. 1: Diffusivity calculated as in reference [H for L = 16 
and several AT, proving the consistency of assumptions also 
for small size (and high gradients) systems. 



in the neighborhood are up and half are down (see e.g. 

aai). 

KS move: consider a diagonal with two opposite spins, 
and exchange them whenever energy is preserved (see 
0). 

The second-neighbours exchange in KS is essential for 
the dynamization of otherwise frozen configurations near 
the cold border, ensuring an effective transitivity in the 
configuration space. Then the evolution rule may be de- 
fined as follows: 

KQ Dynamics: a single KQ step is a sequence oflxi 
randomly alternated Q2R and KS moves on randomly 
chosen sites and diagonals. Such a step defines the nat- 
ural time unit r. 

Besides tests already performed in [5(, the reliability 
of the KQ dynamics has been successfully checked by 
looking at the robustness of the results with respect to 
various perturbations. A meaningful test, for instance, 
consists in a neat change of the randomness criterion in 
the choice of sites and diagonals to be moved. By using 
a RW path (which could be also a physically reasonable 
procedure) we obtained indeed the same results, possibly 
apart the time scale. In all cases, a steady state is easily 
established. 

Another important aspect we have verified is that even 
for small systems (L — 16) with large temperature differ- 
ences (Ti = 0, T2 = 7) the energy flow can be described 
by means of a Fourier-like equation with an energy de- 
pendent diffusivity. Therefore, data reported in Fig. [1] 
should be seen as an improvement of those in the figure 
10 of [5|. This confirms the correctness of the ansatz and 
the reliability of the results presented there also very far 
from local equilibrium, i.e. independently of any refer- 
ence to quasi equilibrium local temperature. Indeed, it is 
worth underlining that, in this microcanonical context, 
and especially for small sizes, the local temperature is 
not definite inside the lattice. Therefore, the appropriate 
quantity characterizing local properties is the mean local 
energy. 
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B. Thermodynamic observables 

Typical quantities considered in [5( are the mean en- 
ergy densities of the columns, or (E x ), where x is the 
column label, and averages run for each x on both time 
and column sites. This may be seen as a particular case 
of a general frame. By averaging at every time t along 
the Y direction only, E x = E x (t) is a discrete time se- 
ries; analogously for the squared magnetization M x (t) of 
the x-th column. All typical statistical features of time 
series, first of all time averages and mean square devia- 
tions, may be easily calculated. As usual, well stabilized 
values out of long runs will be considered equivalent to 
the asymptotic ideal values for all practical purposes. 

An interesting point consists in the systematic compar- 
ison between the Ising model with a heat flow (or IMF, 
for brevity) and the closed Ising model (or CIM), i.e. the 
two dimensional toroidal lattice whose energy (a constant 
of motion in this case) will be fixed with suitable criteria. 
Alternatively, one can compare the IMF and the thermal- 
ized Ising model (or TIM) where the flow is zero because 
the borders are fixed at the same temperature. 

More precisely, such comparisons require the following 
steps: 1 - evaluate the mean energy E x of a particular 
column x in the IMF; 2 - fix equal border temperatures 
for the TIM or the total energy of the system for the CIM 
in such a way that the average energy of any column in 
these systems is equal to E x ; 3 - follow the time evolution 
of the systems (IMF, TIM and CIM) in order to obtain 
three sequences of decorraleted values for the different 
observables (e.g. E$ and M|); 4 - compute statistical 
properties of the obtained time series. 

These comparisons aim to stress the influence (if any) 
of the local flow on physical observables with respect to 
different types of thcrmalized systems. 

Of course, an additional check is the comparison be- 
tween TIM and CIM, which should converge to the same 
behaviour for all observables at least when L — > oo. 



C. Geometrical observables 

In order to give evidence to possible correlations be- 
tween heat flow and configurational features, we need a 
different kind of observables. Such observables have al- 
ready been used to study equilibrium states in spin sys- 
tems (precisely Ising systems, with or without long range 
correlations) proving useful in focusing certain peculiar- 
ities of configurations around the critical phase [ill, E3] • 

The precise definition of these quantities requires the 
formalism of configuration and partition spaces, as briefly 
summarized in the Appendix. However, the main idea 
is the following: consider the configuration of a column 
a = a(x, t) as a discrete periodic array of L binary val- 
ues. A probability measure fi is easily defined on the 
array subsets by the normalized number of nodes in each 
subset. This way, an array a (or more precisely the triple 
constituted by a, fj, and the algebra of subsets) becomes 



a particularly simple example of finite probability space. 
An array may be partitioned into homogeneous clusters 
{Ai,A-2,...,A n } of consecutively aligned spins, and this 
collection a = {A\ , A2, ■■ ., A n } may be seen as an element 
of the "partition space" Zl built on the probability space 
a(x). We have established a correspondence $ between 
configurations and finite measurable partitions, or, more 
explicitely 

a(x,t) = (f>(a(x,t)) . 

In this case, the natural order of the cluster sequence 
identifies the partition by the first Y coordinates of each 
cluster. Shannon entropy, conditional entropy, Rohlin 
and Hamming distances between two arbitrary columns 
are therefore well defined functionals (see Appendix) . We 
shall consider in particular the following observables: 

1. the Rohlin distance at a time t between partitions 
a{x) and a(x+l) associated to consecutive columns 
a(x) and a(a; + 1) of the same system, i.e. 

d R (x, t) = d R (a(x, t), a(x + 1, t)) ; 

this distance is a measure of the non similarity be- 
tween adjacent columns, with regard to the cluster 
distributions; 

2. the Rohlin distance between decorrelated columns 
with the same energy (same label x). This is a 
measure of the non similarity between indepen- 
dent columns. Decorrelated configurations can be 
obtained considering either two distinct systems 
evolving independently, or the same system and an 
evolution time At much larger than the decorrela- 
tion time. Therefore, the fistances we consider are 

dii(a.(x, t),f3(x,t)) or dn(a(x,t), a(x,t + At)) 

respectively. 

3. the Hamming distance between two adjacent 
columns, as in item 1, i.e.: 

d H (x, t) = d H (&(x, t), a(x + 1, t)) ; 

this is another and very different measure of non 
similarity, focusing on pointwise differences, inde- 
pendently of the neighborhoods. Moreover, in this 
case, du (x, t) represents the energy between a col- 
umn and the next one. 

The observables defined above, like the previous thermo- 
dynamic quantities, produce discrete time series, admit- 
ting statistical analysis (means, deviations, etc.). 

III. NUMERICAL EXPERIMENTS 

As anticipated, numerical experiments tend to stress 
the influence of the heat flow on significant observables, 
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by comparing IMF and CIM or TIM. Data, in the follow- 
ing, will refer to both time averages and averages over 
multiple experiments. Time averages extend as usual up 
to stable results. Actually, averages run over 10 5 — 10 6 
sampled values, ensuring an excellent stabilization, as if, 
for all practical purposes, the limit t — + oo had been 
reached. 



A. Energy and Magnetization 

The first quantity we shall consider is the energy den- 
sity along the X direction, or E x , x being the label 
of the array a(a;), the configuration of the x-th col- 
umn. For each column the mean in the Y direction is 
always assumed. Consider a system sampled at times 
to,ti,t2, where the starting to occurs after a suitable 
transient (e.g. 50 to 100 times r for L = 16). Moreover, 
in order to have sufficiently decorrelated configurations, 
At = tk—tk-i > 100t. Several At have been tested. The 
resulting time series {E x (tk)} depends also on L and the 
border temperatures (Xi,T2). Then, for every x there is 
a mean energy density (E x ), and a Mean Square Devia- 
tion F = {E x ) — {E x ) 2 (here F stands for fluctuation). 
Such diagrams are plotted in Fig|2]for L = 16, 32, 64 and 
128 at fixed (T U T 2 ). Here T x = 0.01, T 2 = 4, and the 
same in the following, otherwise differently stated. In the 
same figure, at the prescribed energies (E x ), the fluctua- 
tions of the closed system (CIM), are plotted. Since they 
almost coincide for different sizes, only the case L = 16 
is reported, with the error bars. These diagrams show 
that: 

• discrepancies AF between IMF and CIM, defined 
as 

AF = F IMF - F CIM (2) 

(obviously, this definition may be adapted to vari- 
ous cases and observables) , are especially important 
around the critical energy density E c , in a range 
AE = (E c — Si, E c + 82), with Si very small. More- 
over, AF > 0, i.e. fluctuations are always greater 
for IMF. 

• Both the width of AE and the maximal amplitude 
of AF depend on L. Indeed, as L grows, AE de- 
creases and discrepancies AF slowly shrink. The 
way AE decreases seems faster in fact than the 
correlated way the max | AF | vanishes. 

• By comparing data relevant to different sizes, we 
find that, within AF and sufficiently far from F c , 
AF scales like 1/L. As for the very critical point 
F c , our numerical data do not allow any accurate 
prediction about the behaviour of AF, however, 
they suggest that AF decreases slower than 1/L, 
as L grows. Interestingly, this can be read as a 
weak trace of criticality around E c . 
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FIG. 2: Energy fluctuations of columns vs. their mean energy, 
for CIM of size L = 16 (continuous line) and for IMF of 
different sizes (as shown in the legend). The enlarged window 
shows details in the region E c — 81 , E c + 62 . 

We remark that the neighborhood of a certain column 
undergoing a heat flow becomes more and more indistin- 
guishable from an equilibrium neighborhood as L — ► 00. 
Accordingly, it is plausible that the column properties, 
inasmuch as they are related to the state of its neighbor- 
hood, tend to mimic the equilibrium properties in this 
limit. 

In the same spirit, in Fig [3] we can observe, at fixed 
L = 16, the effect of lowering the difference AT = T%—T\ 
for IMF. The convergence of IMF to CIM is again clear, 
starting from E c — 8\ up to E c + <5 2 , where 8 X is very thin 
and 82 is smaller and smaller as L grows. 

Neatly below E c — 8i, or above E c + S2, the coincidence 
between IMF and CIM is quite good for all L and AT. A 
natural question is the reproduction of the same results 
using a TIM instead of a CIM, i.e. a thermalized system 
with equal border temperatures, such to give suitable 
mean energies for comparisons. As a matter of fact, both 
CIM and TIM give indeed qualitatively equivalent results 
with respect to IMF; however, at the observed sizes, they 
do not coincide (see again Fig. [3]). One expects, of course, 
that only for sufficiently high L's a good agreement will 
take place. 

In general, the observed behaviour confirms a fact al- 
ready noticed in Q, i.e. enlarging L is equivalent to 
zooming on a system with a lower AT, so that the ther- 
modynamic limit should give to every column the same 
features of a system in local equilibrium. Clearly, such 
a zooming property is not an absolute equivalence, since 
a finite size TIM cannot reproduce an infinite size IMF. 
The equivalence refers only to the onset of local equilib- 
rium due to the vanishing of the gradient between left 
and right side of each column. Moreover, critical prop- 
erties could disturb the continuity of this process around 
E c . 

Consider now the squared magnetization M 2 , which 
above the critical energy coincides with the mean square 
deviation of M. For a fixed size (here L = 64), in Fig 0] 
we plot the mean values of M 2 vs. energy: the IMF 



FIG. 3: Energy fluctuations of columns vs. their mean en- 
ergy, for different AT in IMF, and for CIM (continuous line). 
The enlarged window shows details putting in evidence the 
intermediate TIM behaviour. For all systems the size is the 
same L — 16. 



FIG. 5: Hamming distance fluctuations of adjacent columns 
vs. the distance self (or longitudinal energy) for several sizes 
of IMF, showing the progressive convergence to CIM in the 
same region previously indicated. 
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FIG. 4: Mean values of squared magnetization M 2 of columns 
vs. mean energy for IMF (stars) and CIM (circles) . The right 
side of the figure is equivalent to the M-fluctuations. 

diagram is neatly above the CIM diagram in the same 
region previously identified by energy fluctuations, from 
E c —8x up to E c +52- Hence, in the same domain, also the 
magnetization fluctuations are larger in the IMF system. 



B. Metric properties 

The energy between a column and the two adjacent 
ones (X direction) should feel, in principle, the asymme- 
try between left and right neighbourhoods. Clearly, as 
remarked in subsection III C[ such a longitudinal energy 
between close columns coincides with their Hamming dis- 
tance du (see Appendix for definitions), giving this met- 
ric concept also a physical interpretation. In Fig. the 
expected difference between IMF and CIM for this quan- 
tity may be easily recognized, once again in the same 



region previously evidenced by energy density and mag- 
netization. 

A quantity directly related to the configurations, more 
precisely to the cluster distributions, is the Rohlin dis- 
tance du (see Appendix), which may be measured with 
various attitudes. For generic partitions, dji(C,rf) is the 
amount of information necessary to distinguish £ from 
?7, i.e. a measure of their non-similarity. Such a non- 
similarity, in our case, can regard both spatially or tem- 
porally distinct cluster distributions. Since this appears 
deeply related to the variability of configurations, d^ is 
a good candidate, in principle, to be an indicator of the 
influence of a gradient on steady states. First of all, we 
consider couples of adjacent columns, so that the longitu- 
dinal energy dn is a meaningful alternative abscissa. The 
mean values and fluctuations of da vs. du are plotted in 
Fig.[5]and[7]respectively, confirming the larger variability 
of IMF system. 

It would be also interesting to understand if it is pos- 
sible to distinguish systems with or without heat flow by 
looking at a single column. To this end we consider the 
sequence of uncorrelated configurations at times t\ . . .tk, 
calling ckfc = a(x,tk) the corresponding partition for the 
x-th column (see Appendix for details). We calculate 
the numerical sequence of distances: dn{ak, afc+i). Such 
a sequence follows the "novelty creation" along an or- 
bit for every examined column, whereas the previous se- 
quence followed the evolution of an isochronous gradi- 
ent of novelty between adjacent columns. In both cases, 
fluctuations give overall estimates of such dynamic or 
isochronous variability. 

In Fig. [5] we observe the behaviour of time averaged 
dft{ak, at+i) for L = 16,32,64 as a function of the en- 
ergy of the corresponding columns. For clearness, we 
have splitted the comparison in two frames, 16-32 and 16- 
64 respectively. Apart the incidental inversion between 
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FIG. 6: Rohlin distance between consecutive rows vs. Ham- 
ming distance between the same rows. Results pertaining to 
IMF and CIM are compared. 
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FIG. 8: Time averaged Rohlin distances of single columns as 
a function of the average energy of the column itself. Mean 
values for L = 16, 32 (left window) and L = 16, 64 (right 
window) are depicted, as shown in the legend. 
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FIG. 7: Rohlin distance fluctuations vs. the Hamming dis- 
tance for IMF and CIM and two sizes (L = 16, 32). 

IMF and CIM at L = 16, at larger L, IMF-distances 
are greater than the corresponding CIM-distances. Once 
again the larger variability of the system presenting heat 
flow is evidenced. 

We note also that the maximum evolves with L: the 
peak grows logarithmically, as expected, while the peak 
abscissa slowly decreases. 

As to fluctuations, results summarized in Fig. [9] are 
extremely similar to those in the previous Fig. [7J 

Three points have to be stressed: 

• the remarkable likeness between Figures [7] and [5] is 
far from being trivial, since partitions are strictly 
correlated in the former case, uncorrelated in the 
latter; 

• fluctuations are only slightly wider in the uncorre- 
lated cases; 

• the region interested by a discrepancy between IMF 
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FIG. 9: Fluctuations of Rohlin distance for single columns 
vs. their mean energy for IMF and CIM and two sizes (L 
In 32). 

and CIM is again the same, possibly with a small 
shift toward low energies for the left bound. 

A further remark is that nothing would be different in 
Fig. [5] using partitions from two independent systems: 
this confirms the complete decorrelation of configurations 
along an orbit within At. 



IV. CONCLUSIONS AND PERSPECTIVES 

All numerical experiments on the Ising cylindrical 
model converge on the fact that the imposed heat flow re- 
veals in a wider amplitude of fluctuations for local observ- 
ables. Recalling the robustness of present results with 
respect to small perturbations of the dynamics, as re- 
marked in (jll A[) , a natural question arises: how much 
does this behaviour depend on the IMF peculiarities? In 
other terms, would an asymmety in the boundary con- 
ditions, as the left-right temperature difference in our 
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model, be automatically translated into amplified fluc- 
tuations, when imposed on a generic lattice system? If 
so, being a general consequence of spatial asymmetry in 
probabilistic processes, this feature would be very weakly 
related to physics. We would argue, on the contrary, that 
the observed behaviours of IMF vs. CIM or TIM are non 
trivially related to real mesoscopic features of a magnetic 
system. 

First of all, a point stressing the physical meaning of 
our experiments is that the influence of heat flow on ob- 
servables appears to be deeply related to the peculiar 
way an Ising rectangular model passes through the crit- 
ical region. The amplification, as remarked, does not 
regard indeed the whole of a steady state, but only a rel- 
evant neighborhood of the magnetic transition. On the 
contrary, for small values of T± and high values of T2, 
observables in the regions close to the borders are prac- 
tically indistinguishable from those in equivalent equi- 
librium states. This last feature may be understood in 
terms of typical configurations: near the cold border, 
there are indeed sparse spots of one or two sites, mak- 
ing the left and right neighborhoods of the observed col- 
umn practically identical. The same happens near the 
hot border, provided that the temperature is sufficiently 
high to establish a uniform disorder, this time because of 
the irreducible fragmentation into thin clusters. Only in 
the intermediate region there is a meaningful difference 
between left and right sides, reflecting the growth and 
subsequent fragmentation of clusters in the X direction. 
Columns are slices of such clusters, with a shape depen- 
dence on x heavily related the properties of the Ising 
system. 

A further indication that an asymmetry in boundary 
conditions is not sufficient to explain the larger variability 
of IMF is provided by a simple study of the paradigmatic 
model of non interacting diffusion, i.e. RW on a lattice. 
Precisely, by imposing different densities of walkers at the 
borders, it is possible to show that, even in the presence 
of a strong density gradient, fluctuations in the system 
remain unchanged. 

Hence, a purely diffusive RW is too poor to reproduce 
the behaviour we have observed in the Ising model, where 
evidently interaction plays a fundamental role. In the 
same way, the very existence of a critical temperature 
(or energy) , which is certainly related to the observed ef- 
fects, is irreproducible by simple RW. In order to clarify 
the subject, local interactions should be introduced in the 
RW model, mimicking the role of the energy dependent 
diffusivity in IMF. This may be done in several ways, and 
studies in this direction are in progress, as well as tests 
on totally different dynamical systems (e.g. asymmet- 
ric sandpiles). All this will be fully reported in another 
paper. 

Finally, we remark that the relevance of a finite (i.e. 
non infinitesimal) thermal gradient, or the consequent 
vanishing of AE and AF in the thermodynamic limit, do 
not imply that the observed effects are physically mean- 
ingless. There are no reasons indeed to consider finite 



size properties as unphysical. A mesoscopic situation (L 
finite) with peculiar non-equilibrium features could be 
equally or even more interesting from a physical point of 
view. 
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APPENDIX A: APPENDIX: CONFIGURATIONS, 
PARTITIONS SPACES AND DISTANCES 

Let M be a graph with L nodes or sites <ij assuming 
values in an alphabet K. A configuration on M is a whole 
set a = {aj}, aj £ K. It is an element of C = C(M), the 
set of all |K| L possible states of the lattice. For instance, 
if M is a discrete array (as in the case of our columns) 
or a square lattice, and K = { — 1, 1}, this description fits 
Ising-like systems. 

A path, is a sequence of "near" sites, and a connected 
cluster is a set of sites with the same value in K which 
are connected by a path. For general graphs, clusters 
are connected but not necessarily simply connected sets. 
Since every site belongs to a single cluster, clusters Ak 
are disjoint subsets of M and (J fe Ak = M. In other 
terms, the clusters collection is a "finite partition" of M, 
whose subsets {Ak] constitute its "atoms" . The partition 
space Z = Z(M) is the set of all finite partitions of M. 
The correspondence $ : C — > Z between a configuration 
a e C and the clusters partition a = (At, ...,Ax) G Z, 
i.e. a = $(a), is "many to one", because configurations 
generated by permutations in K are mapped into the 
same partition. 

For every subset A of M, let n(A) be the normalized 
number of nodes in A. This defines a probability measure 
/i in the algebra M. of subsets of M. 

For standard operations on partitions in Z(M) classi- 
cal textbooks are e.g. [H, [TJ, [l5|, [l6|]. For applications 
in the spirit of our demands, see also [Til . fl2l [T7L fl8| . 
Here we only recall the definitions of Shannon entropy 
and Rohlin distance. 

Let a = (A\, An) be a partition: its Shannon en- 
tropy H(a) is 

N 

H ( a ) = -5>(^)lnM(^) . (Al) 

The Shannon entropy does not depend on the shapes 
of the atoms, but only on their measures. If = 
(i?i, -Bm") is another partition, shapes implicitly in- 
fluence the conditional entropy of a with respect to (3: 

N M , . , 

H(a\0) = - E E n B k ) In ^ . (A2) 
i=i fc=i k > 

Then, the Rohlin distance da between partitions is de- 
fined by 

d R (a,0)=H(a\0)+H(0\a) . (A3) 
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This makes Z(M.) a metric space. The Rohlin distance 
expresses how different two partitions are. 

If K itself is a metric space (e.g. a numerical set with 
the usual distance between numbers), one can also con- 
sider in C(M) the Hamming distance dn which, for con- 
figurations a and b, is defined by the functional 

d H ( a ,b)=J2\b j -a j | (A4) 

3 

(possibly normalized by dividing by L). In our case, as 
noticed in Section fill Bl the Hamming distance between 
adjacent columns is the energy between them. 

In general, Hamming and Rohlin distances are not di- 
rectly comparable. The former is between configurations, 
and it is sensitive only to actual values of corresponding 
nodes, not to their distribution or neighborhood, whereas 
the latter is between partitions, and therefore it is sen- 
sitive to the cluster shapes. In principle, du and dn 
may give very different information. With a binary al- 
phabet, for instance, complementary configurations have 
maximal Hamming distance (dn — L), while the corre- 
sponding partitions coincide (dn = 0). 

If a configuration a £ C has discrete evolution 

a, Ssl, S a, S a, ... 

one can speak of "configurations orbit" . The correspond- 



ing dynamics S on Z is defined by 

Sa = S $(a) = $ (5a) (A5) 

so that to a configurations orbit there corresponds a par- 
titions orbit {S n a} = {$(5"a)}. Clearly, the probability 
measure fi in Z is not preserved by S, because clusters 
do not evolve in themselves but are redefined at every 
step by the pointwise dynamics in C. However, we are 
not interested here in indicators requiring a preserved 
measure, such as Kolmogorov-Sinai entropy or Lyapunov 
exponents. 

Real valued observables F or F, in C(M) or Z(M), 
give rise to "time series" {fk} = {F(S k a)} or {f k } — 
{F(S k a)}. Such time series are typical objects of our 
investigations. 

This formalism applies in principle to every kind of 
lattices and discrete dynamics. Note however that when 
M is a one dimensional array, as in the case considered 
here, the Rohlin distance is essentially simpler than in the 
two-dimensional case, because of the geometrical nature 
of the atoms contours: points in the former ca se, p ossi- 
bly cumbersome paths in the latter (see e.g. [TTI. 1 1 2t fTsj] ) . 
For the Hamming distance, on the contrary, the compu- 
tational complexity would be almost the same. 
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